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We consider the problem of designing an optimal quantum detector that distinguishes unambigu- 
ously between a collection of mixed quantum states. Using arguments of duality in vector space 
optimization, we derive necessary and sufficient conditions for an optimal measurement that maxi- 
mizes the probability of correct detection. We show that the previous optimal measurements that 
were derived for certain special cases satisfy these optimality conditions. We then consider state 
sets with strong symmetry properties, and show that the optimal measurement operators for dis- 
tinguishing between these states share the same symmetries, and can be computed very efficiently 
by solving a reduced size semidefinite program. 

PACS numbers: 03.67.Hk 



I. INTRODUCTION 

The problem of detecting information stored in the 
state of a quantum system is a fundamental problem in 
quantum information theory. Several approaches have 
emerged to distinguishing between a collection of non- 
orthogonal quantum states. In one approach, a mea- 
surement is designed to maximize theprobability of cor- 
rect detection HSHSHHISlESIi. A more 
recent approach, referred to as unambiguous detection 
[ISIIlOEQHEIIlllIlisto design a measure- 
ment that with a certain probability returns an incon- 
clusive result, but such that if the measurement returns 
an answer, then the answer is correct with probability 
1. An interesting alternative approach for distinguishing 
between a collection of quantum states, which is a com- 
bination of the previous two approaches, is to allow for 
a certain probability of an inconclusive result, and then 
maximize the probability of correct detection [Tj| [jjj, HOI • 

We consider a quantum state ensemble consisting of m 
density operators {pi, 1 < i < m} on an n-dimensional 
complex Hilbert space Tl, with prior probabilities {pi > 
0, 1 < i < m}. A pure-state ensemble is one in which 
each density operator pi is a rank-one projector \4>i)((f>i\, 
where the vectors \4>i), though evidently normalized to 
unit length, are not necessarily orthogonal. Our problem 
is to design a quantum detector to distinguish unambigu- 
ously between the states {pi}. 

Chefles [l5| showed that a necessary and sufficient con- 
dition for the existence of unambiguous measurements 
for distinguishing between a collection of pure quantum 
states is that the states are linearly independent. Nec- 
essary and sufficient conditions on the optimal measure- 
ment minimizing the probability of an inconclusive result 
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for pure states were derived in ^t|- The optimal mea- 
surement when distinguishing between a broad class of 
symmetric pure-state sets was also considered in |l7|. 

The problem of unambiguous detection between mixed 
state ensembles has received considerably less attention. 
Rudolph et al. 21] showed that unambiguous detection 
between mixed quantum states is possible as long as one 
of the density operators in the ensemble has a non-zero 
overlap with the intersection of the kernels of the other 
density operators. They then consider the problem of un- 
ambiguous detection between two mixed quantum states, 
and derive upper and lower bounds on the probability of 
an inconclusive result. They also develop a closed form 
solution for the optimal measurement in the case in which 
both states have kernels of dimension 1. 

In this paper we develop a general framework for un- 
ambiguous state discrimination between a collection of 
mixed quantum states, which can be applied to any 
number of states with arbitrary prior probabilities. For 
our measurement we consider general positive operator- 
valued measures 0, , consisting of m + 1 measurement 
operators. We derive a set of necessary and sufficient 
conditions for an optimal measurement that minimizes 
the probability of an inconclusive result, by exploiting 
principles of duality theory in vector space optimization. 
We then show that the previous optimal measurements 
that were derived for certain special cases satisfy these 
optimality conditions. 

Next, we consider geometrically uniform (GU) and 
compound GU state sets 0, 0, H3| , which are state sets 
with strong symmetry properties. We show that the op- 
timal measurement operators for unambiguous discrim- 
ination between such state sets are also GU and CGU 
respectively, with generators that can be computed very 
efficiently by solving a reduced size semidefinite program. 

The paper is organized as follows. In Section [HI we 
provide a statement of our problem. In Section lTTTl we de- 
velop the necessary and sufficient conditions for optimal- 
ity using Lagrange duality theory. Some special cases are 
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considered in Section HV1 In Section we consider the 
problem of distinguishing between a collection of states 
with a broad class of symmetry properties. 

II. PROBLEM FORMULATION 

Assume that a quantum channel is prepared in a quan- 
tum state drawn from a collection of mixed states, rep- 
resented by density operators {pi, 1 < i < to} on an n- 
dimensional complex Hilbert space 7i. We assume with- 
out loss of generalit y th at the eigenvectors of pi , 1 < i < 
77i, collectively span 31] Ti. 

To detect the state of the system a measurement is 
constructed comprising m + 1 measurement operators 
{IL,, < i < to} that satisfy 

LTi > 0, < i < to; 

ET=o^ = i- (i) 

The measurement operators are constructed so that ei- 
ther the state is correctly detected, or the measurement 
returns an inconclusive result. Thus, each of the opera- 
tors Hi, 1 < i < to correspond to detection of the corre- 
sponding states pi, 1 < i < to, and IIo corresponds to an 
inconclusive result. 

Given that the state of the system is pj, the proba- 
bility of obtaining outcome i is Tr(pjIF). Therefore, to 
ensure that each state is either correctly detected or an 
inconclusive result is obtained, we must have 

Ti(pjUi) = rjiSij, l<i,j<m, (2) 

for some < rji < 1. Since from IIo = I — Y7L\ H») 
(0 implies that Tr(pillo) = 1 — rji, so that given that the 
state of the system is pi, the state is correctly detected 
with probability rji , and an inconclusive result is returned 
with probability 1 — r )j. 

It was shown in ^H] that for pure-state ensembles con- 
sisting of rank-one density operators pi — \4>i)(4>i\, (0 
can be satisfied if and only if the vectors \<pi) are linearly 
independent. For mixed states, it was shown in [2l| that 
(0 can be satisfied if and only if one of the density op- 
erators pi has a non-zero overlap with the intersection of 
the kernels of the other density operators. Specifically, 
denote by ICi the null space of pi and let 

Si = nf^j^JCj (3) 

denote the intersection of JCj,l < j < m,j ^ i. Then 
to satisfy the eigenvectors of 11^ must be contained 
in Si and must not be entirely contained in ICi. This 
implies that ICi must not be entirely contained in Si. 
Some examples of mixed states for which unambiguous 
detection is possible are given in plf . 

If the state pi is prepared with prior probability pi, 
then the total probability of correctly detecting the state 
is 

TO 

P D = J^Pi^iPi^i)- ( 4 ) 

i=l 



Our problem therefore is to choose the measurement op- 
erators IIj,0 < i < to to maximize Pp, subject to the 
constraints and 

Tr( ft n i ) = 0, l<i,j<m,i^j. (5) 

To satisfy IL; must lie in Si defined by ©, so that 

Ili = PilliPi, 1 < i < m, (6) 

where Pi is the orthogonal projection onto Si. Denoting 
by Qi an n x r matrix whose columns form an arbitrary 
orthonormal basis for Si, where r = dim(5i), we can 
express P, as Pi — 0^0*. From © and Q we then have 
that 

Hi = OiAiQ*, l<i<m, (7) 
where Aj = 0*11^0^ is an r x r matrix satisfying 
A l > 0, 1 < i < to; 
I7=i 0^,6* < /. (8) 
Therefore, our problem reduces to maximizing 

m 

Pd = ^PiTtip&i&iQX), (9) 

i=l 

subject to (JSJ. 

To show that the problem of © and © does not de- 
pend on the choice of orthonormal basis Oi, we note that 
any orthonormal basis for Si can be expressed as the 
columns of ^i, where = 0;t/i for some r x r unitary 
matrix C/j. Substituting ^i instead of Qi in @ and JSJl, 
our problem becomes that of maximizing 

m m 

P D = ^ftTrO^Ai**) = ^^Tr^e.A^G*), (10) 

i=l i=l 

where A- = UiAiU*, subject to 

A l > 0, 1 < i < to; 

Since A l > if and only if A- > 0, the problem of (fTU|) 
and Ullfl is equivalent to that of © and JSJl. 

Equipped with the standard operations of addition and 
multiplication by real numbers, the space B of all Hermi- 
tian 77. x 7i matrices is an n 2 -dimensional real vector space. 
As noted in [2l| , by choosing an appropriate basis for B, 
the problem of maximizing Pd subject to © can be put 
in the form of a standard semidefinite programming prob- 
lem, which is a convex optimization problem; for a de- 
tailed treatment of semidefinite programming problems 
see, e.g., 0> EH EH 113 • By exploiting the many well 
known algorithms for solving semidefinite programs |27| , 
e.g., interior point methods |32]| 0)E^|j the optimal mea- 
surement can be computed very efficiently in polynomial 
time within any desired accuracy. 

Using elements of duality theory in vector space op- 
timization, in the next section we derive necessary and 
sufficient conditions on the measurement operators Ilj = 
QiAiO* to maximize P D of © subject to ©■ 
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III. CONDITIONS FOR OPTIMALITY 
A. Dual Problem Formulation 

To derive necessary and sufficient conditions for opti- 
maiity on the matrices Aj we first derive a dual problem, 
using Lagrange duality theory psj . 

Denote by A the set of all ordered sets II = 
{IT, = e.AiO*}™! satisfying © and define J(II) = 
Y^iLi PiTr(pi9iAi0*). Then our problem is 



max J(U). 

neA v ; 



(12) 



We refer to this problem as the primal problem, and to 
any II £ A as a primal feasible point. The optimal value 
of J(II) is denoted by J. 

To develop the dual problem associated with (|12l) we 
first compute the Lagrange dual function, which is given 
by 

9(Z) = 



mm < - } ^Tr(p t 0A9*) 



Tr z hTe^e*-/ 



= mm|^Tr(A i e*(^-PiPi)e i )-Tr(Z)|, (13) 

where Z > is the Lagrange dual variable. Since 
Aj > 0, 1 < i < to, we have that Tr(A 2 A) > for 
any X > 0. If A is not positive semidefinite, then we 
can always choose A, such that Tr(AiA) is unbounded 
below. Therefore, 

= f -Tr(Z), A, > 0, 1 < i < to, Z > 0; (u 
y ^ ' y — oo, otherwise, ^ ' 

where 

A i = Q* i (Z-p iPi )Q i , l<i<m. (15) 
It follows that the dual problem associated with 11211 is 
nunTr(Z) (16) 

subject to 

Qt( z -Pi/H)&i>0, l<i<m; 



Z>0. 



(17) 



Denoting by L the set of all Hermitian operators Z such 
that e*{Z - p l p i )Q i > 0, 1 < i < to and Z > 0, and 
defining T(Z) = Tr(Z), the dual problem can be written 
as 



minT(Z). 

zer 



(18) 



We refer to any Z £ L as a dual feasible point. The 
optimal value of T(Z) is denoted by T. 



B. Optimality Conditions 

We can immediately verify that both the primal and 
the dual problem are strictly feasible. Therefore, their 
optimal values are attainable and the duality gap is zero 
i.e., 



In addition, for any II 

Zer, 

T(Z) - J(U) = 



J = T. (19) 
= {IT = QiAi@*i}il i G A and 



\t=i 

> 0, 



UnZ 



(20) 



where Tl = I - YT=i ©i^iO* > 0. Note, that gQJ 
can be used to develop an upper bound on the optimal 
probability of correct detection J. Indeed, since for any 
ZeT, T(Z) > J(II), we have that J < T(Z) for any 
dual feasible^. 

Now, let Ui = 9 4 A 4 e*,l < i < m and n = / - 
Si=i denote the optimal measurement operators that 
maximize © subject to |[5}, and let Z denote the optimal 
Z that minimizes (|16|) subject to l|17|) . From l|19|l and 
(|20|l we conclude that 



Tr \^U&*{Z -piPijQi + fk 



Z =0. 



(21) 



Since A; > 0, Z > and e*(Z-p l/ 9 i )9. i > 0, 1 < i < to, 
is satisfied if and only if 



ZU = 

e* i (z- PiPi )e i A i = o, 



1< i < m. 



(22) 
(23) 



Once we find the optimal Z that minimizes the dual 
problem (|16fl . the constraints Il22[l and 123|) arc nec- 
essary and sufficient conditions on the optimal mea- 
surement operators IT. We have already seen that 
these conditions are necessary. To show that they 
are sufficient, we note that if a set of feasible mea- 
surement operators IL satisfies (|22|l and (|23|l . then 
Tr(YZi^i @ i( Z -PiPi) @ i+^oZ) = so that from 

ED, J(n) = T(Z) = J. 

We summarize our results in the following theorem: 

Theorem 1. Let {pt,l < i < m} denote a set of density 
operators with prior probabilities {pi > 0, 1 < i < m}, 
and let {8i,l < i < m} denote a set of matrices such 
that the columns of O^ form an orthonormal basis for 
Si = HijLi j^lCj , where JCi the null space of pi. Let 
A denote the set of all ordered sets of Hermitian mea- 
surement operators II = {IL}™ that satisfy 11; > 0, 
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$^,-_oHj — 1 '> an d Tr(pjTli) = 0,1 < i < m,i ^ j and 
let r denote the set of Hermitian matrices Z such that 
Z > 0, Q^(Z— PiPi)Qi, 1 < i < m. Consider the problem 
maxngA </(n) and the dual problem mmzgr T(Z), where 
J(J1) = YZiPiTHPi^i) and T(Z) = Tr(Z). Then 

1. For any Z ET and U e A, T(Z) > J(II). 

2. There is an optimal II, denoted II, such that J = 
J(fi) > J(n) for any U £ A; 

3. There is an optimal Z , denoted Z and such that 
f = T{Z) < T(Z) for any Z E T; 

4-f= J; 

5. Necessary and sufficient conditions on the optimal 
measurement operators 11, are that there exists a 
Z e T such that 

Zlto = (24) 
e*(Z - Pip)QiKi = 0, 1 < i < m, (25) 

where Hi = OiAi®*, 1 < i < ra, and Aj > 0. 

6. Given Z , necessary and sufficient conditions on the 
optimal measurement operators Ilj are 

zn = (26) 
Q*{Z - PiPi)QiAi = 0, 1 < i < m. (27) 

Although the necessary and sufficient conditions of 
Theorem ^ are hard to solve, they can be used to ver- 
ify a solution and to gain some insight into the optimal 
measurement operators. In the next section we show that 
the previous optimal measurements that were derived in 
the literature for certain special cases satisfy these opti- 
mality conditions. 



IV. SPECIAL CASES 

We now consider two special cases that where ad- 
dressed in |2l|, for which a closed form solution exists. In 
Section II V Al we consider the case in which the spaces Si 
defined by © are orthogonal, and in Section llV Bl we con- 
sider the problem of distinguishing unambiguously be- 
tween two density operators with dim(Si) = 1, 1 < i < 2. 



A. Orthogonal Null Spaces Si 

The first case we consider is the case in which the null 
spaces Si are orthogonal, so that 



where Pi is an orthogonal projection onto Si. It was 
shown in |2l| that in this case the optimal measurement 
operators are 



u l = p l = e t e*, i < i < m. 



(29) 



In Appendix lAl we show that the optimal solution of the 
dual problem can be expressed as 



z = 22piPiPiPi. 



(30) 



It can easily be shown that Z and IT of l|30|) and l|29|) 

satisfy the optimality conditions of Theorem ^ 



B. Null Spaces of Dimension 1 

We now consider the case of distinguishing between 
two density operators p± and p2, where S± and 52 both 
have dimension equal to 1. In this case, as we show in 
Appendix [51 the optimal dual solution is 

(diPu d 2 -d 1 \f\ 2 <0; 

Z=ld 2 P 2 , di-rf 2 |/| 2 <0; 
[d 2 {@2 + sQ 2 L )(Q 2 + sQ^)*, otherwise, 

(31) 

where Pi is an orthogonal projection onto Si, Q 2 is a 
unit norm vector in the span of Oi and & 2 such that 
6^ = 0, and 



d i =p i e* i p i Q i) l<z<2; 
/ = e|9 i; 

e=(e 2 L ye 1 . 



(32) 



The optimal measurement operators for this case were 
developed in |2l| . and can be written as 



where 



n 1 = p 1 ,n 2 = o, d 2 -di|/| 2 <0; 
n 1 = o,n 2 = p 2 , di-rf 2 |/| 2 <o ; 

III = «i-Pi,n 2 = a 2 P 2 , otherwise, 

(33) 



d 2 \tl- 



a 2 



i-l/l 2 



dl|/| 2 



(34) 



PiPj = Sij, 1 < i,j,< m, 



(28) 



We now show that Z and II of l|31|l and (|33|) satisfy 
the optimality conditions of Theorem ^ To this end we 
note that from (|33l) . 

'A 1 = 1,A 2 = 0, d 2 -di|/| 2 <0; 
(A i }? =1 =^Ai=0,Aa = l, rfi-rf 2 |/| 2 <0; 
Ai = ai, A 2 = a 2 , otherwise. 

(35) 
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From (EU-GUl we have that if d 2 - di|/| 2 < 0, then 
&t(Z - pipi)e x Ai = di - ejpiPiGi = 0; 

e^(z-p 2 p 2 )e 2 A 2 -o; 

Zn = Z(J - III) = diGiej - dxQi&t = 0. (36) 
Similarly, if di - d 2 \f\ 2 < 0, then 

01(^-^)0x^1 = 0; 

95 - P2P2 )e 2 A 2 = d 2 - e^ P2 p 2 e 2 = o ; 
z% = z(i - fi 2 ) = d 2 e 2 ej - d 2 e 2 e^ = o. (37) 

Finally, if neither of the conditions d\ — rf 2 |/| 2 < 0, d% — 
^i|/| 2 < hold, then 

Qt(Z - p lPl )Q 1 A 1 = 



(d 2 (r + e*s)(r + e*s)*-d 1 )- 



d 2 \f\' 



1-1/1^ 



d 2 \f\ 2 




l-l/l 2 



0. 



(38) 



e^(z-p 2/92 )e 2 A 2 = (e* 2 ze 2 -d 2 ) 



i-l/l 5 



o, 



(39) 



and 



zrio = z - zni - zn 2 

= z- AiZGiej - A 2 ze 2 e 2 . (40) 

To show that ZIIo = 0, we note that 

ze^i = d 2 (\f\ 2 + s *ef*)e 2 e* 2 
+ d 2 ( s \f\ 2 + ss *er)eie* 2 

+ d 2 (e*/ + s*|e| 2 )0 2 2 L * 

+ d 2 (se*f + ss*\e\ 2 )QiQi*, (41) 



and 



ze 2 e 2 = d 2 e 2 e 2 + d 2S e 2 L e 2 . 



(42) 



Substituting JIT} and (|4*2|) into iflUjl. and after some al- 
gebraic manipulations, we have that 



OPTIMAL DETECTION OF SYMMETRIC 
STATES 



We now consider the case in which the quantum state 
ensemble has symmetry properties referred to as geomet- 
ric uniformity (GU) and compound geometric uniformity 
(CGU). These symmetry properties are quite general, 
and include many cases of practical interest. 

Under a variety of different optimality criteria the op- 
timal measurement for distinguishing between GU and 
CGU state sets was shown to be GU and CGU respec- 
tively 0)IS0j0|- I n particular it was shown in that 
the optimal measurement for unambiguous detection be- 
tween linearly independent GU and CGU pure-states is 
GU and CGU respectively. We now generalize this result 
to unambiguous detection of mixed GU and CGU state 
sets. 



VI. GU STATE SETS 

A GU state set is defined as a set of density operators 
1 < i < m} such that pi = UipU* where p is an 
arbitrary generating operator and the matrices {Ui, 1 < 
i < to} are unitary and form an abelian group Q [a. l29|~. 
For concreteness, we assume that U\ = I. The group 
Q is the generating group of S. For consistency with 
the symmetry of S, we will assume equiprobable prior 
probabilities on S. 

As we now show, the optimal measurement operators 
that maximize the probability of correct detection when 
distinguishing unambiguously between the density oper- 
ators of a GU state set are also GU with the same gen- 
erating group. The corresponding generator can be com- 
puted very efficiently in polynomial time. 

Suppose that the optimal measurement operators that 
maximize 



j({n i }) = gir( ft n i ) 



(44) 



subject to JSJ and (JSJ arc U;, and let J = J({Ui}) = 
Y^iLi Tr(pilF). Let r(j, i) be the mapping from TxX to X 
with J = {1, . . . , to}, defined by r(j, i) = k if U*U t = Uk- 

Then the measurement operators IT 



-0) 



n 



0) 



UjU r (j t i)Uj and 



ZU = Z- AiZOiej - A 2 Z9 2 e 2 = O. (43) > 0,1 < i < to and 



/ — Y^aL x n,- for any 1 < j < to are also optimal. 

Indeed, since Ilj > 0, 1 < i < m and Y^iiLi < I, 
U) 



Combining l|3tj|) - (|43[l we conclude that the optimal mea- 
surement operators of |2lj | satisfy the optimality condi- 
tions of Theorem 2] 



CO _ 



u* < U 3 U* = I. 



(45) 
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Using the fact that pi = UipU* for some generator p, 

m 

i=l 
m 



m 

i=l 
J. 



(46) 



Finally, for I 7^ i, 
Tr (^fif 



TV (UipU^U^ A U*) 
Tr iu s pU*U r{j ^) 



Tr (jO s II fc 
0, 



(47) 



where U s = U*Ui and U k = U*Ui and the last equality 
follows from the fact that since I 7^ i, s 7^ k. 

It was shown in 0, 0] that if the measurement 

operators IF^ are optimal for any j, then {IF = 

(l/m)YT=i n| j) , 1 < i < m} and TT = / - E^IF subject to the constraints 



VII. CGU STATE SETS 

A CGU state set is defined as a set density operators 
{pik, l<«<^,l<fc<?"} such that p ik = Ui<f)kU* for 
some generating density operators {p k , 1 < k < r}, where 
the matrices {Ui.1 < i < Z} are unitary and form an 
abelian group Q ja.l23|. A CGU state set is in general not 
GU. However, for every k, the operators {pit, 1 < i < i} 
are GU with generating group Q. 

Using arguments similar to hose of Section IVTl and 18] 
we can show that the optimal measurement operators 
corresponding to a CGU state set can always be chosen 
to be GU with the same generating group Q as the orig- 
inal state set. Thus, to find the optimal measurement 
operators all we need is to find the optimal generators 
IF. The remaining operators are obtained by applying 
the group Q to each of the generators life. 

Since the optimal measurement operators satisfy IUfc = 
t^IW,l < i < 1,1 < k < r and p lk = U lPk U*, 
Ti(p ik U ik ) = Tr(p k H k ), so that the problem reduces 
to the maximization problem 



llfcfco ■ — ' 



(50) 



fc=i 



are also optimal. Furthermore, IF = UiUU? where 

We therefore conclude that the optimal measurement 
operators can always be chosen to be GU with the same 
generating group Q as the original state set. Thus, to 
find the optimal measurement operators all we need is to 
find the optimal generator II. The remaining operators 
are obtained by applying the group Q to II. 

Since the optimal measurement operators satisfy IF = 
UiUU*, 1 < i < m and Pl = U^U*, TrO^IF) = Tr(pII), 
so that the problem @ reduces to the maximization 
problem 



IFc > 0, 1 < k < r; 

E-=iELi^n fc c/* fe </; 

Tr(n fcPlfe ) = 0, l<A<r,2<t<l. 



(51) 



Since this problem is a (convex) semidefinite program 
ming problem, the optimal generators IF can be com- 
puted very efficiently in polynomial time within any de- 
sired accuracy |24L I26L l27| . Note that the problem of E~ 



max Tr (oil), 



(48) 



where B is the set of n x n Hcrmitian operators, subject 
to the constraints 



II > 0; 
YZ 1 UiUU* < I; 
Tr(np,) = 0, 2 < i < to. 



(49) 



The problem of (|48|l and (|49|l is a (convex) semidefinite 
programming problem, and therefore the optimal II can 
be computed very efficiently in polynomial time within 
any desired accuracy [2J, |2(| |23 , for example using the 
FMI toolbox on Matlab. Note that the problem of 
and l|49l) has n 2 real unknowns and m + 1 constraints, 
in contrast with the original maximization problem @ 
subject to © and (JjJ which has ran 2 real unknowns and 



m 



1 constraints. 



and l|51|) has rn 2 real unknowns and lr+1 constraints, in 
contrast with the original maximization which has Irn 2 
real unknowns and (lr) 2 + 1 constraints. 



VIII. CONCLUSION 

We considered the problem of distinguishing unam- 
biguously between a collection of mixed quantum states. 
Using elements of duality theory in vector space opti- 
mization, we derived a set of necessary and sufficient 
conditions on the optimal measurement operators. We 
then considered some special cases for which closed form 
solutions are known, and showed that they satisfy our op- 
timality conditions. We also showed that in the case in 
which the states to be distinguished have strong symme- 
try properties, the optimal measurement operators have 
the same symmetries, and can be determined efficiently 
by solving a semidefinite programming problem. 

An interesting future direction to pursue is to use the 
optimality conditions we developed in this paper to de- 
rive closed form solutions for other special cases. 
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APPENDIX A: PROOF OF (1551) 



APPENDIX B: PROOF OF STDi 



To develop the optimal dual solution in the case of 
orthogonal null spaces, let = [0i 02 ■■■ 0m], and 
define a matrix ©^ such that [0 -1 ] is a square, uni- 
tary matrix, i.e., [0 0^] [0 0^] = /. Denoting 
Z = [0 Q 1 -] Y [0 0^]*, the dual problem can be ex- 
pressed as 



minTr (J© 0^] Y [0 Q 1 -]' 



(Al) 



subject to 



0* [0 0^] Y [0 0^]*©, > 0|ftPi6i, 1 < i < m; 

Y > 0. (A2) 

Using the orthogonality properties of ©i and 0^, the 
problem of (|A1() and l|A2|l is equivalent to 



subject to 



min Tr(Y") 

Y 



Yi > QiPiPiQi, 1 < i < m: 
Y > 0, 



(A3) 



(A4) 



where 



Y 



r. 



F 2 



y 

1 n, 







(A5) 



Since Tr(F) = E™i Tr ( r 0' a solution to (O subject 
to EH is 



y 



>2 



y 



(A6) 



where 



y = ®*p l p l Q l , l<i<m. 



(A7) 



Then, 



Z = [0 x ]y[0 ± ]* = ^p i P i p i P i , (A8) 



To develop the optimal dual solution Z for one- 
dimensional null spaces, we note that Z lies in the space 
spanned by 0i and ©2- Denoting by a matrix whose 
columns represent an orthonormal basis for this space, Z 
can be written as Z = ©y©*, where the 2x2 matrix Y 
is the solution to 



subject to 



min Tr(y ) 

Y 



$;y$i > di; 
$*y$ 2 > d 2 ; 

y > o. 



(Bl) 

(B2) 
(B3) 
(B4) 

we 



Here $ 4 = @*0j and di — Pi<d*pi<di for 1 < i < 2. 

To develop a solution to l|BT|) subject to (|B2|) - (|B4 
form the Lagrangian 



C = Tr(y) - £) 7i (*?y$ i - - Tr(Xy), (B5) 

i=i 

where from the Karush-Kuhn- Tucker (KKT) conditions 
[3fj| we must have that 7, > 0, X > 0, and 

7i ($?Y$ i -d i ) =0, £ = 1,2; (B6) 
Tr(Ay) = 0. (B7) 

Differentiating £ with respect to Y and equating to zero, 



2 

£ 

i=l 



7i $i$J - X = 0. 



(B8) 



If A = 0, then we must have that I = X) i=l7 i$»$*, 
which is possible only if $1 and $2 are orthogonal. 
Therefore, X 7^ 0, which implies from 1|B7[) that 1|B4|) 
is active. Now, suppose that only (|B4|I is active. In this 
case our problem reduces to minimizing Tr(j/*y) whose 
optimal solution is y — 0, which does not satisfy l|B2|) 
and CEU). 

We conclude that at the optimal solution 1|B4|) and at 
least one of the constraints i|B2(l and (|B3(I are active. 
Thus, to determine the optimal solution we need to de- 
termine the solutions under each of the 3 possibilities: 
only l|B2|) is active, only (|B3(1 is active, both (|B2(1 and 
(|B3|) are active, and then choose the solution with the 
smallest objective. 

Consider first the case in which (|B2(1 and l|B4|) are ac- 
tive. In this case, Y = yy* for some vector y, and without 
loss of generality we can assume that 



$*y = di. 

To satisfy (|B9|) . y must have the form 



(B9) 



EDJ. 



(B10) 



8 



where is a unit norm vector orthogonal to $1, so that 
$*<J)J- — 0, and s is chosen to minimize Tr(Y). Since, 



Tr(y) = y*y — di + \s\ 2 , 



(Bll) 



s = 0. Thus, F = and Tr(F) = di. This 

solution is valid only if l|B3|l is satisfied, i.e., only if 

$* 2 Y<f> 2 = d 1 \f\ 2 >d 2 . (B12) 

Here we used the fact that 

= 9^69*9! = 9^9i = /, (B13) 

since 99* is an orthogonal projection onto the space 
spanned by 9i and 9 2 . 

Next, suppose that (|B3J| and l)B4|l are active. In this 
case, Y = yy* where without loss of generality we can 
choose y such that 



and 



y 



d 2 $ 2 + S$o 



(B14) 



(B15) 



where $> 2 1S a UIU t norm vector orthogonal to $2, and 
s is chosen to minimize Tr(Y). Since Tr(Y) = d 2 + \s\ 2 , 
s = 0, and Tr(Y) = d 2 - This solution is valid only if l|B2|) 
is satisfied, i.e., 



^\Y^ 1 =d 2 \f\ 2 >d 1 . 



(B16) 



Finally, consider the case in which l)B2jl - (|B4|l are ac- 
tive. In this case, we can assume without loss of gener- 
ality that = V^2- Then, 



y = yd 2 <& 



2 • 



where s is chosen such that 

$;y$i = d u 



(B17) 



(B18) 



and Tr(Y) = y*y is minimized. Now, for y given by 
p7)l . 

Y = d 2 $2$ 2 + |s| 2 $2 L $ 2 L * + 

+s-v/^*2 $ 2 + s*\/d2$2$£*, ( B1 9) 



so that 



$*Y$! = d 2 |/| 5 



l*| 2 |e| 2 + 



^Z + ^el 2 , 



d 2 se*f + ^/d 2 s*f*e 
(B20) 



where we defined 9^ = 9$^! an d e an d / are given by 
(|32[) . Therefore, to satisfy l|B18[) . s must be of the form 



for some <£>. The problem of IjBljl then becomes 
1 

mm — — 

<p \e\ z 

which is equivalent to 



max3?{e^/} 



Since 

R{e*V} < |e^/| =|/|, 
the optimal choice of <p is = /*/|/|, and 



e* 



^7 



- 1 



(B21) 



(B22) 



(B23) 



(B24) 



(B25) 



For this choice of s, 



Tr(F) = d 2 + \s\ 



= d 2 1 + 



A 

= a. 



I/I 2 / 



- 1 



(B26) 



Clearly, a > d 2 . Therefore, to complete the proof of 
(|31|l we need to show that a > d\ . Now, 



|e| 2 (a-di) 



|e| 2 (d 2 - d x ) + |/f 



fdl 



l/l 



= (1 - |e| 2 H + (|e| 2 + |/| 2 )d 2 - 2v^v^|/| 
= (\f\Vd~1-Vd~2) 2 

> 0, (B27) 

where we used the fact that 

|e| 2 + |/| 2 = 9t9 2 9^9 1 + 9t9 2 L (9 2 L )*9 1 

= e* 1 Q 1 = 1, (B28) 

since 9 2 9*. + 9^(9^)* is an orthogonal projection onto 
the space spanned by 9i and 9 2 . 
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